Quantum interference in coherent molecular conductance 
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Coherent electronic transport through individual molecules is crucially sensitive to quantum inter- 
ference. Using exact diagonalization techniques, we investigate the zero-bias and zero-temperature 
conductance through 7r-conjugated annulene molecules (modeled by the Pariser-Parr-Pople and 
Hubbard Hamiltonians) weakly coupled to two leads. We analyze the conductance for different 
source-drain configurations, finding an important reduction for certain transmission channels and 
for particular geometries as a consequence of destructive quantum interference between states with 
definite momenta. When translational symmetry is broken by an external perturbation we find an 
abrupt increase of the conductance through those channels. Previous studies concentrated on the 
effect at the Fermi energy, where this effect is very small. By analysing the effect of symmetry 
breaking on the main transmission channels we find a much larger response thus leading to the 
possibility of a larger switching of the conductance through single molecules. 



The possibility of achieving controlled quantum trans- 
port through single molecules has become a reality 
as seen from various successful attempts in different 
systems 0, 0, Q. The electronic transport through sin- 
gle 7r-conjugated molecules has been studied in several 
theoretical 0, 0, 0, 0, S] and experimental^, [13, El El 
works. This important step leads not only to miniatur- 
ization of electronic devices but also to the possibility 
of taking advantage of new physical properties in these 
systems. 

A fundamental physical property in quantum trans- 
port is quantum interference, which could be a handle to 
control conductance through such systems. Constructive 
and destructive interference play a crucial role which, for 
most cases, is per se independent of the structure and 
composition of the molecular bridges to the leads. 

Even though changing from constructive to destructive 
interference might be possible for certain geometries and 
molecules, disrupting the destructive interference by per- 
turbing the molecule seems to be more dramatic. This 
has been proposed in two recent theoretical works @, 0] 
for the so-called Quantum Interference Effect Transis- 
tors (QuIET) based on single annulene molecules, includ- 
ing benzene. In the first paper, the systems are mod- 
eled by a Pariser-Parr-Pople (PPP) Hamiltonian 13[ us- 
ing the Ohno parametrization (l3 | and solved using the 
self-consistent Hartree-Fock (HF) approximation, while 
the second work resorts to ab-initio calculations within 
the LDA and ab-initio approximations. The conductance 
in both cases is then calculated for the strong coupling 
limit using the non-equilibrium Green's function |l5| and 
the Landauer-Biittiker formalism 16] . Both studies con- 
centrate on the equilibrium conductance at zero bias and 
gate voltage (V g — 0, i.e. the off-resonance tunnelling 
regime) which coincides with the Fermi energy of the 



leads. However, annulenes have a gap at this energy due 
mainly to the strong Coulomb interactions present in the 
molecule (N = An + 2 annulenes like benzene, are closed- 
shell molecules with no level at zero energy even in the 
absence of interactions). The conductance at zero gate 
voltage is finite, albeit small, only for strong coupling to 
the leads[l3]- For weak coupling the conductance will 
be zero and will be appreciable only through the main 
channels of the molecule which are a few eV away from 
the Fermi energy. By analysing the QuIET at the main 
transmittance channels the switching effect will be much 
more pronounced and robust as we will show below. 

In this work we analyze the linear conductance of a 
series of annulenes at finite gate voltages and several 
source-drain lead configurations. For leads set up at 
180 degrees (the A configuration in Fig. 1), the sys- 
tem presents completely constructive interference in all 
transmission channels. However, for other lead config- 
urations, we find different behaviours depending on N. 
For N = An + 2 annulenes we find a reduced conduc- 
tance through particular channels due to a nearly de- 
structive interferences explained below. When the molec- 
ular translational invariance is broken by means of elastic 
scattering or decoherence factors, we find an abrupt in- 
crease in the conductance through these channels due to 
disruption of destructive interference. This increase is 
much larger than was found in previous works since the 
focus is set on a resonance transmission channel through 
the molecule. 

However, for N = An annulenes the effect is even more 
striking: shifting the drain lead by one lattice site from 
the A to the B configuration (see Fig. 1) reduces the 
transmission to zero. This is, in fact due completely to 
the destructive interference for N = An molecules as ex- 
plained below. A translational symmetry breaking per- 
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FIG. 1: Schematic representation of the molecules studied nu- 
merically in the presence of a gate voltage V g , small hybridiza- 
tion to the leads t' and two source-drain lead geometries. 



turbation will have, as in this case, a huge effect, since 
the conductance through these channels will increase dra- 
matically from zero to an appreciable value. Hence, such 
molecules could have important applications. Previous 
theoretical studies analyze the effect of electron correla- 
tion on the structure stability of carbon rings exhibiting 
competing many-body effects of aromaticity, Jahn- Teller 
distortion and dimerization[l8l Hij. These studies indi- 
cate that, under certain conditions, some of the systems 
considered below could be stabilized. In addition, ring 
geometries of quantum dots are also potential devices 
where such a behaviour could be observed. 

In this work we consider an annular 7r-conjugated 
molecules with N sites, weakly connected to non- 
interacting leads in the A or B configurations (Fig. [I}. 
We consider the following Hamiltonian for the system 
consisting of two leads weakly coupled to the molecule: 
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The first term describes the isolated 7r-conjugated 
molecule, modeled by the PPP Hamiltonian . with on- 
site energy given by a gate voltage V g : 

firing = ~t C L c 3o +Y U i ( n H ~ \ ) ( n *T ~ \ 
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where the operators c\ a (ci a ) create (annihilate) an elec- 
tron of spin a in the tt orbital of the Carbon atom at site 
i, n are the corresponding number operators and (• • • ) 
stands for bonded pairs of Carbon atoms. The intersite 
interaction potential Vij is parametrized so as to inter- 



polate between U and e /r^ in the limit r\ 
In the Ohno interpolation, Vij is given by 
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where the distance 



is m 



A. 



The standard Hubbard 
parameter for sp 2 Carbon is Ui — 11.26 eV and hopping 
parameter t for r = 1.397 A is 2.4 eV [20( , and all energies 
are in eV. The second term in Eq.(l) corresponds to two 
tight-binding semi-infinite chains for the left and right 



leads. The third term in Eq. (JTJ) describes the coupling 
of the edge sites of the left and right leads with sites L 
(left) and R (right) of the system respectively. When 
the ground state of the isolated ring is non-degenerate, 
and the coupling t' between the leads and the ring is 
weak (t' <C t), equilibrium conductance at K can be 
expressed to second order in t' in terms of the retarded 
Green's function for the isolated ring between sites i and 
j: Gf'jfa). For an incident particle with energy ui = 
-It cos k and momentum ±fc, the transmittance reads 
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where e 1(a) H = t' 2 G R LL{RR) (u;), i» = f J G» fl ( U ), 
play the role of a correction to the on-site energy at 
the extremes of the leads and an effective hopping be- 
tween them respectively. The conductance is G = 
(2e 2 /h)T{fi, V g ), where \x is the Fermi level, which we set 
to zero (half-filled leads) and V g enters implicitly through 
the Green's functions. 

In our study, the interacting system with N sites is 
solved exactly using numerical techniques (Lanczos or 
Davidson) to diagonalize huge matrices, and the Green's 
functions are obtained straightforwardly The correla- 
tions are, thus, treated in an exact manner. 

We now analyze the two possible scenarios for annu- 
lenes with N = An + 2 and N = An with two terminals. 
For the first case, the allowed total momentum quan- 
tum numbers are k — 2rit / {An + 2) = rix j{2n + 1), with 
r an integer, while for the second case the allowed mo- 
menta are k = Irix j An — m/2n. For a two-terminal set 
up, wave functions travelling through both branches of 
the molecule will interfere producing different interfer- 
ence patterns depending on the positions of the leads. 
The phase difference will be momentum times the differ- 
ence in the lenghts of the two trajectories (in units of the 
C-C separation): A0 = kAx. 

Specializing to benzene (N = 6), for leads in the 
"para" position, Ax = and the waves are in phase, in- 
terfering constructively (Fig. [2]). However, in the "meta" 
position (Ax = 2), the interference will depend on the 
k value of the particular channel. For the highest oc- 
cupied molecular orbital (HOMO) and lowest unoccu- 
pied molecular orbital (LUMO) the phase differences will 
be A(f> = 2ir/3 and Atfi — Att/3 respectively and the 
interference will reduce the amplitude in these chan- 
nels. For benzene configuration of the leads produc- 
ing a completely destructive scenario for these channels 
does not exist. Previous experimental and theoretical 
works [23|, [23| noted that the transmission through ben- 
zene in the "meta" configuration is much lower than that 
through the "para" position. 

There is no transmission channel at V g = since the 
density of states of the bare molecule at that energy is 
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FIG. 2: Top: Transmittance through a benzene molecule 
for the "para" (A substitution in Fig. 1) (broken line) and 
"meta" (B substitution in Fig. 1) (full line) configurations 
as a function of the gate voltage measured from the Fermi 
energy. Here t' = 0.4. Bottom: Same for the "meta" config- 
uration in the presence of different on-site potentials which 
break the translational invariance. Here a finite Lorentzian 
width r\ — 0.03 and t' — 0.1 have been taken for visualization 
purposes. Inset: evolution of the area under the transmission 
peak through the HOMO as a function of local energy. 

zero. A detailed analysis of quantum interference effects 
in benzene was carried out for the large coupling regime 
[25| where interesting effects appear at zero energy and 
finite bias [26|, (2?| at weak coupling. 

What would happen if the translational symmetry is 
broken by an external perturbation? This question was 
first addressed in [6| by introducing a local energy (£), 
at one site in benzene, the real part of E would produce 
elastic scattering and its imaginary part, decoherence. 
In that study, the focus was at the Fermi energy of the 
leads (set to zero) where the observed effect is small. 
However, by studying the effect of external perturbations 
on the main transmittance channels such as the HOMO 
and LUMO, we find a much larger response. We show 
these results in the bottom part of figure 2, where an 
additional diagonal energy is added to the site to the right 
of the B ( "meta" ) position (the effect is not qualitatively 
dependent on this position). It is clearly seen that, in 
this case, the small peak corresponding to the HOMO 
level develops and grows as the local energy is increased, 
disrupting the translational symmetry responsible for the 
destructive interference (see inset). 

For larger An + 2 rings, the momentum of the HOMO 
level approaches 7r/2 and the interference becomes more 
destructive. This is shown for a 10 site annulene in figure 
3, where the peaks for the B position of the leads are 
much smaller than those corresponding to benzene (Fig. 
2, top). Growth of the HOMO peak with a local potential 
energy is larger in this case than in benzene. 

However, for annulencs with 4n atoms, depending on 
the total momentum of the particular channel, the inter- 



FIG. 3: Same as figure 2, for N = 10 annulene. 

ference can be totally destructive. For these annulenes 
the interference will be completely constructive in the A 
configuration and totally destructive if one lead is shifted 
by one site (B configuration, see Fig. 1). In these cases 
we can observe the emergence of a transmission channel, 
if the translational invariance is disrupted. 

The effect of symmetry breaking by a local external 
perturbation is much more striking in annulenes with 
An sites. As discussed above, we expect a new peak to 
emerge when the translational symmetry is broken. In 
figure 4 we show results for conductance through an 8- 
site annulene represented by the PPP model. In the top 
figure, the HOMO and LUMO channels, fully formed for 
the A configuration, completely vanish in the B case due 
to quantum interference. When a local perturbation is 




FIG. 4: Top: Transmittance through an [8]-annulene molecule 
for the A (broken line) and B (full line) configurations as a 
function of the gate voltage measured from the Fermi energy. 
Here t' = 0.4. Bottom: Transmittance for the B configuration 
in the presence of a local perturbation of the hopping of the 
form i(l — ip) which breaks the translational invariance (n = 
0.03 and t' — 0.1). Inset: evolution of the HOMO integral as 
a function of ip. 
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FIG. 5: Top: HOMO integral of a [12]-annulene molecule for 
the A (broken line) and B (full line) configurations as a func- 
tion of the locally perturbed hopping t(l — <p) (left panel) and 
the dimerization parameter S (right panel). Bottom: Trans- 
mittance for the B configuration in the presence of dimeriza- 
tion. 



configuration (Fig. 1) the interference will be nearly de- 
structive. The larger the molecule the closer the phase 
difference between different paths is to ir and hence 
the interference becomes more destructive. For [An]- 
annulcncs. on the other hand, the interference for these 
configurations can be completely destructive. We have 
shown the effect of symmetry-breaking perturbations on 
different transmission channels. For [4ro]-annulenes the 
effect is striking since the transmission changes from zero 
to a large finite value for small perturbations. These ef- 
fects should be seen in substituted annulenes and should 
also be appreciable in transport through rings of quan- 
tum dots. 

We thank A. A. Aligia and B. Alascio for useful discus- 
sions. This work was done in the framework of projects 
PIP 5254 of CONICET, PICT 2006/483 of the ANPCyT 
and the ARG/RPO-041/2006 Indo- Argentine collabora- 
tion. SR thanks DST for support through JC Bose fel- 
lowship. 



applied to the molecule, in this case a different hopping 
on the bond to the right of the B lead, t(l — <p), the ab- 
sent peak emerges as shown in the inset of the bottom 
figure. A similar effect occurs also for a local potential 

Finally, in figure 5 we present results for a [12]- 
annulcnc. In this case also the HOMO and LUMO peaks 
for the B configuration are non-existent due to destruc- 
tive interference. Here we show how the HOMO peak 
arises as a function of a different hopping parameter in 
one lead of the form t(l — ip) (top left) and as a function 
of dimerization, 6, for transfer integral for t(l — (— 1) J <5) 
of the bond between sites i and i + 1 (top right). The 
transmittance in the A configuration doesn't change in 
the presence of dimerization, while in case B, the HOMO 
and LUMO channels develop suddenly already for a small 
S (bottom panel) . Dimerization of the annulenes lifts the 
two-fold degeneracy of the HOMOs as well as the LUMOs 
in large annulenes. Thus, the condition for destructive 
interference between the two paths is lifted and one ob- 
serves transmission through these channels. 

In conclusion, by analysing the resonant conductance 
through the HOMO and LUMO channels in the weak 
lead-molecule coupling regime we find a strong depen- 
dence on the source-drain configuration and on the 
molecular geometry due to quantum interference. This 
effect is more robust and striking than in the strong- 
coupling case since in the latter the destructive interfer- 
ence occurs only at zero gate voltage and can be masked 
in an experiment, while in the former a whole channel can 
appear or disappear depending on the geometry. We'd 
also like to stress that these results are not affected by 
molecular vibrations at room temperature since modes 
that can cause decoherence are excited at temperatures 
higher than 500^0. 

In [An + 2]-annulenes like benzene in the B ("meta") 
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